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ABSTRACT 


The  purpose  of  this  paper  is  to  study  some  recent  applications  of 

the  n by  dn  LCP  solvable  by  a parametric  principal  pivoting  algorithm 

(PPP  algorithm).  First,  it  is  shown  that  by  analyzing  the  n by  dn  LCP 

we  could  study  the  problem  of  solving  a system  of  equations  and  the 

(nonlinear)  complementarity  problem  when  the  function  involved  is  separable. 

Next,  we  examine  conditions  under  which  the  PPP  algorithm  is  applicable  to 

a general  LCP,  and  then  present  examples  of  LCP's  arising  from  various 

applications  satisfying  the  conditions ; included  among  them  is  the  n by  dn 

LCP  with  a certain  p-property.  Finally  we  study  a special  class  of 

' 

n by  dn  LCP's  which  do  not  possess  the  P-property  but  to  which  the  PPP 

' 

algorithm  is  still  applicable;  a major  application  of  this  class  of 
problems  is  a certain  economic  spatial  equilibrium  model  with  piecewise 
linear  prices. 

\ 


1.  INTRODUCTION 


...  zSt.  . 


Given  a positive  integer  d,  an  n-vector  q,  n by  n matrices  M*,...,Md 
1 d-1 

and  positive  n-vectors  f , the  n by  dn  linear  complementarity 

1 d 

problem  (n  by  dn  LCP)  is  that  of  finding  n-vectors  w,x  , . . . ,x  satisfying 
the  conditions : 

(1.1a)  w=q+lJ=1Mixi 

(1.1b)  w > 0,  0 < x1  < f1  , i=1 , . . . ,d-1 , xd  > 0 

T 1 i i T i+1 

(1.1c)  w X = 0,  (f1  - X1) 1 X1  =0,  i=1 , . . . ,d-1  . 

When  d=1 , (1.1)  reduces  to  the  (ordinary)  LCP,  which  is  to  find  m-vectors 
w and  x such  that 

(1.2)  w = r + Lx,  w > 0,  x > 0,  wTx  = 0, 

where  r and  L are  given  m-vector  and  m by  m matrix,  respectively.  On  die 
other  hand,  (1.1)  can  be  converted  into  an  ordinary  LCP  (1.2)  with  a dn  by  dn 
L (see  [12J). 

The  n by  dn  LCP  was  first  studied  by  one  of  the  authors  in  [12]  who 
characterized  the  existence  and  uniqueness  of  the  solution  of  the  problem  by 
a certain  P- property  (see  Section  2).  It  was  also  shown  that  when  the 
problem  has  this  property,  then  the  solution  can  be  obtained  by  applying  a 
parametric  principal  pivoting  algorithm  (abbreviated  as  PPP  algorithm  in  the 
following:  see  Section  3)  to  the  LCP  converted  from  the  original  n by  dn  LCP. 
A noteworthy  point  here  is  that  the  matrix  L in  this  LCP  is  neither  a P-matrix 


2 


nor  positive  semi -definite  but  still  the  PPP  algorithm  is  shown  to  be  applicable. 
As  noted  in  [12]  applications  leading  to  an  n by  dn  LCP  with  the  p-property 
include  a nonlinear  analysis  of  certain  reinforced  concrete  and  other  structures 
(see  Kaneko  [13]),  and  a strictly  convex  quadratic  program  with  bounded 
variables  which  by  itself  has  many  applications  (see  [6,  25]). 

The  purposes  of  this  paper  are  (i)  to  examine  further  the  applicability 
of  the  PPP  algorithm  to  the  n by  dn  LCP,  and  (ii)  to.  study  some  recent 
applications  of  the  n by  dn  LCP,  solvable  by  the  PPP  algorithm. 

In  the  next  section,  we  shall  show  that  two  closely  related  problems, 
i.e.,  (i)  the  problem  of  solving  a system  of  piecewise  linear  equations 
defined  on  a cubic  subdivision  and  (ii)  the  piecewise  linear  complementarity 
problem  defined  on  a cubic  subdivision,  are  equivalent  to  the  n by  dn  LCP. 

This  implies,  in  particular,  that  one  could  analyze  these  problems  by 
considering  the  corresponding  n by  dn  LCP. 

After  discussing  the  applicability  of  the  PPP  algorithm  in  some  general 
terms  in  Section  3,  we  shall  study  in  the  fourth  and  final  section  another 
application  of  the  n by  dn  LCP  arising  from  the  computation  of  equilibrium 
prices  in  an  economic  spatial  model.  In  this  model,  a certain  commodity  is 
transshipped  among  various  regions  where  its  prices  are  given  as  piecewise 
linear  functions  of  the  local  net  imports.  (Of  course,  the  piecewise 
linearity  could  arise  as  an  approximation  of  more  general  nonlinearities.) 

One  wants  to  determine  a set  of  prices  and  commodity  flows  such  that  certain 
spatial  equilibrium  conditions  are  satisfied.  In  most  practical  cases,  the 
n by  dn  LCP  obtained  in  this  model  does  not  possess  the  P-property; 
however,  we  shall  show  that  the  problem  can  still  be  solved  by  the  PPP 
algorithm. 

The  special  case  of  the  equilibrium  model  where  the  prices  are  linear  has 
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been  studied  by  many  researchers  (see  [1,  10,  28,  31]).  In  particular,  one 
of  the  authors  has  shown  in  [28]  that  the  PPP  algorithm  can  be  applied  to 
solve  the  linear  model  and  that  the  computational  procedure  is  greatly 
simplified  by  taking  advantage  of  some  special  structures  of  the  problem. 

The  linear  model  has  a number  of  applications  (see  Polito  [30],  Takayama  and 
Judge  [32]  and  references  therein)  but  suffers  the  obvious  limitation  because 
of  the  linear  assumption  on  the  prices.  Presumably,  the  piecewise  linear 

model  offers  a broader  scope  of  applicability. 

To  our  knowledge,  there  has  been  no  explicit  study  of  the  piecewise 
linear  model.  The  general  model  (where  prices  are  general  nonlinear  functions) 
has  been  treated  in  [ 29  ] where  necessary  and  sufficient  conditions  for  equilibrium 
are  derived  in  terms  of  a generalized  complementarity  problem  over  polyhedral 
cones,  and  in  [20]  where  a fixed  point  algorithm  is  suggested  to  compute 
the  equilibrium  prices  under  some  general  setting. 


2.  THE  n BY  dn  LCP  AND  RELATED  PROBLEMS 

In  this  section  we  shall  study  some  relationships  of  the  n by  dn  LCP 

(1.1)  with  the  following  two  problems: 


(2.1)  Find  x € Rn  such  that  h = G(x),  where  h€Rn,  and  G:  Rn  — > Rn  is 
continuous  and  piecewise  linear  on  a cubic  subdivision  of  Rn  ; 


(2.2)  Find  w € Rn  and  x€Rn  such  that 


w = s + F(x),  w > 0,  x > 0,  w x = 0, 

where  s €Rn,  and  F:  R™  — > Rn  is  continuous  and  piecewise  linear 
on  a cubic  subdivision  of  R^  . 

After  giving  a definition  of  a cubic  subdivision  of  Rn  (or  R^) , we  shall 
show  that  the  three  problems,  (1.1),  (2.1)  and  (2.2)  are  equivalent  in  the 
sense  that  any  one  of  them  can  be  transformed  into  either  of  the  other  two. 

The  requirement  that  a function  be  piecewise  linear  on  a cubic  subdivision 
is  highly  restrictive;  for  instance,  a (nonlinear)  function  must  virtually 
be  separable  in  order  that  it  can  be  approximated  by  a piecewise  linear 
function  on  a cubic  subdivision.  Thus,  the  scope  of  application  of  the 
problems  (2.1)  and  (2.2),  is  admittedly  severely  limited.  However,  investigation 
of  these  problems  seems  to  be  justified  because  (i)  in  many  practical  situations 
nonlinear  functions  representing  the  problem  are  in  fact  separable  (see  e.g., 
problems  described  in  [9,  21,  34,  35]  in  addition  to  well-known  engineering 
problems) , and  (ii)  the  simple  structure  of  the  problem  leads  to  sharp 
theoretical  results  and  more  efficient  computation  of  the  solution.  Some 
results  can  be  found  in  recent  papers  by  Kojima  [16]  and  Todd  [33]  on  the 
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computation  of  a solution  of  the  system  of  equations  (2.1)  using  fixed  point 


algorithms . 


Let  positive  integers  . . . .d^  scalars  fj  , j=1,...,n,  and  positive 
scalars  f^  , i=1,...,d  - 1,  j=t,...,n  be  given.  Let  T be  the  set  of  all  y 
of  the  form  y = (Y.  .....  y)  where  y.  € {0, 1 , . . . ,d  } , j=1,...,n.  For  each 

J J 

y 6 T,  define 

Yj-1  _ Yj 

(2.3)  C(y)  = { x € Rn : £ f!  < r < I f!  if  y.  >1,  and  x.  < f°  if  Y,  = 0} 

i=0  J J i-o  J J J ~ J J 


where  by  convention  we  set  f^J  = + « for  each  j.  We  shall  call  the  partition 
(C(y)  : y 6 T]  of  Rn  a cubic  subdivision  of  Rn. 

A cubic  subdivision  of  R*\  the  nonnegative  orthant  of  Rn,  can  be  similarly 

T 

defined.  Let  d^,...,dn  be  positive  integers  and  f*  positive  scalars, 
i=1 , . . . ,d^  - 1 , j=1 , . . . ,n.  Let  T+  be  the  set  of  y = (y^ , . . . .y^)  with 

y.  €{l,...,d  ) , j=1,...,n.  For  each  y € f+  define 


(2.4) 


V • Yj 

D(y)  = { x € R.“  : £ f J < x.  < £ f J , J-l , . . . ,n  } , 

,*  -n  J J “ J 


where  by  convention  f®  = 0 and  f ® for  each  j. 

A function  G : R > R (or  F : r”  > Rn)  is  called  piecewise  linear 
on  a cubic  subdivision  if  it  is  affine  on  each  cell  C(y)  (or  D(y), 
respectively)  of  the  subdivision.  For  the  sake  of  simplicity,  we  shall  assume— ^ 
in  the  following  that  d = d^  for  all  j.  We  shall  denote  by  fi  the  n-vector 


1 There  is  no  loss  of  generality  in  doing  so;  in  general,  let  d be  the 

largest  d^  and  for  k with  d > d^,  one  may  add  d - d^  dummy  subdivision(s) 

with  respect  to  the  k-th  coordinate. 
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whose  j-th  component  is  , i=0,...,d. 

Before  proceeding  further,  we  would  like  to  observe  the  fact  that  the 

complementarity  conditions  (along  with  nonnegativity)  in  the  problem  (1.1) 

provide  an  equivalent  statement  of  the  condition  x€D(y)  in  the  problem  (2.2). 

To  be  more  specific,  let  [d(y):  y ^ ^ 3 be  a cubic  subdivision  of  Rn.  If 

V • y . - 1 + . 

x€D(Y),  then  by  setting  x = f i =1 , . . . ,Y . -1 , x.J  = x.  - £.jn  f*  and 

J J J J J i=0  j 


xj  = o , i = for  each  j=1 , . . . ,nj  it  follows  that  {x\...,xd}  , 

i i i T 

where  x = (x^ . . . ,x^)  , i=1,...,d,  satisfies  the  conditions  (1.1b)  - (1.1c) 

T 1 Id 

(except  for  w >0  and  w x =0).  Conversely,  if  {x  ,...,x  } satisfies 

(1.1b)  - (1.1c),  then  x = x^+  . . . +xd  € D(y) , where  for  each  j = 1 , . . . ,n, . y^  - 1 

is  the  largest  index  i such  that  x*  * f*  . Similarly,  if  {c(y)  : Y € F}  is 

a cubic  subdivision  of  Rn,  then  z€c(y)  is  characterized  by  the  following 
0 d 

conditions:  z - z + ...  + z and 


z°  < f°  , 0 < zl  < fl  , i=1 , . . . ,d-1 , zd  > 0 , 

(2.5) 

(f1  - z1)1  z1+1  = 0,  i = 0,1,...,d-1  . 

2.1  Proof  of  Equivalence  of  Problems  (1.1),  (2,1)  and  (2.2) 

First  of  all,  we  shall  show  the  equivalence  between  (1.1)  and  (2.2). 

To  prove  that  the  n by  dn  LCP  (1.1)  can  be  written  in  the  form  (2.2), 

it  suffices  to  show  that  q + £d_^  x^  is  a continuous  function  of 

x = x^+  ...  +xd  which  is  affine  on  D(y)  for  each  y in  T+  . This  is  obvious 

2/ 

because  we  have  — 

2 For  a matrix  K,  K.^  denotes  its  j-th  column. 


r 


*1 


I; 

I it 

I 


ilk 


-r 


.i  i 


i=1 


(2.6) 


d 

n 

= ( 

q + 

E 

E . x 

cr 

if 

•i-i  • 

i=1 

j=>  J 

j 

j-1 

E 

M*. 

f1 

Yi 

+ M.  (x . 

V 

- E 

i=1 

J 

J 

j 3 

i=1 

■ -1 

J 

E 

(M* 

Yi  i 

M.-bfb  + 

[m!J, 

i=1 

j 

J 3 

n d 


A i 


j=1  1=1  'j  j 


Y_ 

. . ,M.n]x 

n 


To  show  the  other  implication,  Let  the  piecewise  linear  complementarity 


problem  (2.2)  be  given  and  let  {D(y)  : Y € T+}  be  the  associated  cubic  subdivision 


„n 


of  R+  . For  each  Y^F+  , let 
(2.7)  F (x)  = a (v)  + A(y)x, 


for  some  a (Y)  € Rn  and  A(Y)€RnXn  for  each  x6  D(y).  Further  let 
M*  = A(yS,  where  y*  = for  i=1,...,d.  Then,  it  is  not  difficult 


to  show  (see  [14])  that  for  each  y in  , we  have  that 


Y Y 

(2.8)  A(y)  « [M.!  , . . . ,M.n  ] . 

i n 


Also  it  follows  from  some  elementary  calculations  (see  [14])  that 


(2.9) 


n Yj“^  . Y. 

a(Y)  = b + E E [M^.  - M.J  ]f,  , 

3 3 3 


j=1  i=1 


1, 


where  b = a(y  ).  By  (2.8),  (2.9),  and  by  reversing  the  argument 
in  (2.6),  we  can  conclude  that  (2.2)  is  converted  to  the  n by  dn  LCP  (1.1) 
with  q = s +b. 


It  is  well-known  that  a complementarity  problem  car.  be  transformed 
into  the  problem  of  solving  a system  of  equations;  i.e.,  the  transformation 


* 


from  (2.2)  into  (2.1)  is  known.  However,  the  transformation  must  be  such 
that  the  resulting  equations  are  piecewise  linear  on  a cubic  subdivision 
of  R . This  can  be  achieved  by  using  the  device  found  in  Eaves  [ 7 ] and 
Megiddo  and  Kojima  [23].  To  be  specific,  let  the  problem  (2.2)  be  given 
and  let  [D(y)  : y€F  } be  the  associated  cubic  s>  bdivision  of  Rn  . Then 
the  problem  is  transformed  into  (2.1)  with  h = - s and 


(2.10)  G(z)  = F(P(z))  - P(z)  + z , z €Rn  , 


where  P(«)  is  the  projection  of  Rn  onto  R^  , i.e.,  P(z)  is  the  "plus  part"  of  z. 
It  is  not  difficult  to  see  that  G is  piecewise  linear  on  a cubic  subdivision 

{C(y)  : y€T},  with  C(y)  as  defined  in  (2.3),  where  for  j=1 , . . . ,n,f*?  = 0 

and  fj  is  the  same  as  that  defining  £d(y)  : y € for  i=1,...,d-1. 

Finally,  to  complete  the  proof  of  equivalence  of  the  three  problems, 
we  show  that  (2.1)  can  be  transformed  into  the  form  (1.1).  Let  the  problem 
(2.1)  be  given  and  let  the  associated  cubic  subdivision  of  Rn  be 
(c(y)  : y€T}.  Further,  let  the  function  G be  expressed  as 


(2.11)  G(z)  = a(y)+A(y)z 


for  each  z€C(y).  Then,  using  the  same  argument  as  before  showing  that  (2.2) 
is  transformed  into  (1.1),  the  problem  (2.1)  can  be  rewritten  as 


(2.12) 


i i 

h = b + I N z 


z°  < f°  , 0 < z1  < f1  , i=1 , . . . ,d-1 , zd  > 0 


,ei  i.T  i+1  n n , , 

(f  -z)  z =0,i=0,...  ,d-1 , 
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where  b = a(y°)  and  N1  = A(yX)  with  yL=  (i  , . . . ,i)T  , i=0,...,d.  We  need 

,0  d 

to  assume  that  N (or  N ) is  nonsingular.  Then  by  letting 


,0  0 i i . , 

w=f  - z ,x  = z , i=1,...,d 


(2.13) 


q = (N°)_1(b-h)  + fU  and  M1  = (NU)~ 'n1  , i=1,...,d, 


we  see  that  (2.12)  is  transformed  into  the  n by  dn  LCP  (1.1).- 


3/ 


2.2  Existence  and  Uniqueness  of  Solution 

Let  M*  € RnXn  , 1=1,..., d.  The  n by  dn  matrix  (M* , . . . ,Md)  is  said  to  have 
the  P-property  if  for  each  y€r+,  the  matrix  A(y)  as  defined  in  (2.8)  is  a 
P-matrix  (i.e.,  all  of  its  principal  minors  are  positive).  The  following 
was  proved  in  [12], 


Theorem  1 . Given  f € R°  , f1  > 0,  i=1 , . . . ,d-1 , the  n by  dn  LCP  (1.1)  has  a 
unique  solution  for  every  q6Rn  if  and  only  if  (M1,...,Md)  has  the  P-property. 

Recently  some  results  have  been  obtained  on  the  existence  and  uniqueness 
of  the  solution  of  a system  of  piecewise  linear  equations  (see  e.g.,  [17, 

18,  19,  22]).  When  the  underlying  subdivision  is  cubic,  then  the  following 
necessary  and  sufficient  condition  holds. 


Theorem  2 (Kojlma  C 1 7l ) . The  system  of  piecewise  linear  equations  (2.1)  on 
a cubic  subdivision  has  a unique  solution  for  every  h?Rn  if  and  only  if  the 

Jacobians  of  G on  all  the  cells  in  the  subdivision  are  nonzero  and  of  the 
same  sign. 


3.  Eaves  and  Lemke  [8  ] have  recently  announced  a scheme  by  which  one  can 
transform  a system  of  piecewise  linear  equations  into  a LCP;  it  is  not 
clear  to  the  authors  of  this  paper  whether  or  not  their  scheme  produces 
the  same  LCP  (or  n by  dn  LCP)  as  that  resulting  from  our  procedure 
described  above. 


It  turns  out,  as  one  might  expect,  that  these  two  results  are  equivalent  in 
the  sense  that  they  imply  each  other;  we  shall  prove  this  equivalence  in 
the  remainder  of  this  section. 

Let  the  n by  dn  LCP  (1.1)  be  given  and  let  [d(y)  : y€r+)  be  the 
associated  cubic  subdivision  of  . Transform  (1.1)  into  the  problem 
(2.1)  via  the  problem  (2.2)  as  we  described  in  the  previous  subsection. 
Namely,  first  convert  (1.1)  into  (2.2);  the  function  F in  the  resulting 
problem  is  of  the  form  F(x)  = a(y)  + A(y)x  for  x€D(y),  where  A(y)  is  as 
given  in  (2.8).  The  problem  (2.1)  is  then  obtained  by  using  the  function  G 
in  (2.10):  this  G is  piecewise  linear  on  a cubic  subdivision  {C(y)  : y G F}  . 
To  prove  that  Theorem  2 implies  Theorem  1,  we  need  to  show  that  the  Jacobians 
of  the  G on  all  the  cells  have  the  same,  nonzero  sign  if  and  only  if  the 
matrix  (M^ , . . . ,Md)  in  (1.1)  has  the  P-property. 

Now,  for  each  y€r,  let  y+  be  the  element  of  T+  which  is  obtained  from  y 
by  replacing  each  zero  component  by  unity.  To  be  more  specific,  let  3 be 
the  set  of  j with  y^  = 0 and  6 its  complement  in  fl,...,n}  . Thus 

y+  = (y|,  . ,.,y*)  satisfies:  y"!"  = 1 for  j€f3  and  yl-  = y^  for  j € 6 . Then, 

4/ 

the  Jacobian  matrix  of  G on  C(y)  is  written  as  — (after  an  appropriate 
principal  rearrangement): 


(2.14) 


A(y+)e 


My+), 


where  I is  the  n by  n identity  matrix.  Clearly,  the  determinar":  of  this 
matrix  is  the  same  as  that  of  A(y+)g  ^ if  6 0 and  unity  if  6 = 0 . This 

implies  that  the  set  of  Jacobians  of  G consists  of  the  determinant  of  A(y+) 
and  all  of  its  principal  minors  for  every  y € F+  , plus  unity.  The  condition 

4.  For  a matrix  K and  index  sets  f and  Tj,  KP  „ denotes  the  submatrix  of  K 
with  rows  and  columns  indexed  by  § and  T|  respectively. 


i 


in  Theorem  2 is  then  equivalent  to  that  A(y+)  is  a p-matrix  for  each  y+  in 
T , or  that  (M\...,Md)  has  the  p-property. 

T 

To  show  that  Theorem  1 implies  Theorem  2,  let  the  problem  (2.1)  with 
a cubic  subdivision  {C(y)  : y € r}  of  Rn  be  given.  Let  G(z)  be  expressed 
as  in  (2.11).  To  apply  Theorem  1,  we  need  to  transform  (2.1)  into  an 
n by  dn  LCP.  Recall  that  this  process  required  that  the  matrix  A(y°) 
with  y°  = (0, . . . ,0)  be  nonsingular.  To  prove  Theorem  2,  however,  we 
can  assume  that  all  A(y)  are  nonsingular,  "for  free".  This  is  because 
it  follows  from  a well-known  result  (see  e.g.,  [19])  that  if  some  A(y) 
is  singular,  then  the  problem  (2.1)  can  not  have  a unique  solution 
for  all  h € Rn  . 

We  shall  consider  only  the  case  where  det  A(y^)  > 0;  the  other  case  is 
similar.  Let  N1  = A(yi)  with  y1  = (i,...,i)  for  i = 0,...,d  and  let 
M = (N^)  , i = 1,...,d.  Then,  m\  . . . ,Md  define  the  n by  dn  LCP  obtained 

from  (2.1).  The  condition  in  Theorem  2 is 


(2.15) 


det  A(y)  >0  all  y € T , 


and  we  need  to  show  that  this  is  equivalent  to  that  (m\  . . . ,Md)  has  the 
P-property. 

To  do  so,  for  any  given  y € I",  let  y+  be  the  element  of  F+  associated 
with  y in  the  sense  that  y+  is  obtained  from  y by  replacing  each  zero 
component  by  unity.  Note  that  if  we  let  3 = [j  : y^  = 0}  and  6 its 
complement,  then 


det  (N^)  ^ A(y)  = det 


<(nV'a<y+))IS6  "j 


0 ((N°) " 1 A(y+))fi 6 j 


det  ((N0)'  1 A(y+)) 


' 


! 


12 

From  this  relationship  and  the  fact  that  det  (N^)  ^ > 0,  it  follows 

that  (2.15)  is  equivalent  to 


(2.16)  For  each  y £T,  det ( (N^)  ^A(y+))gg>  0,  where  6 = {j  : y^  > 1} 


and  y"T  = V . > f 1 or  j € 6 , and  v.  = 1 , for  j £ 6 . 

Now,  it  is  not  difficult  to  show  that  (2.16)  is  indeed  equivalent  to  the 
fact  that  the  matrix  (N^)  V(y+)  is  a P-matrix  for  every  y+€F,  or  that 

(M1 , . . . ,Md)  has  the  P-property. 


f I 
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3.  A PARAMETRIC  PRINCIPAL  PIVOTING  ALGORITHM 


In  this  section  we  shall  examine  the  applicability  of  a parametric 
version  of  Graves'  principal  pivoting  algorithm  (ppp  algorithm;  see 
Graves  [11]  for  the  original  version  and  Cottle  [3  ] for  its  parametric 
version)  to  the  (ordinary)  LCP  and  some  n by  dn  LCP's. 

For  m-vectors  r and  p,  m by  m matrix  L the  parametric  linear 
complementarity  problem  is  that  of  finding  m-vectors  w and  x such  that 


(3.1) 


w = r + Xp  + Lx,  w > 0 , x > 0,  w x = 0 , 


for  each  value  of  scalar  X in  a given  interval.  One  approach  for  solving 
an  ordinary  LCP  (1.2)  is  to  choose  an  appropriate  p and  apply  the  PPP 
algorithm  to  the  parametric  LCP  (3.1);  if  the  algorithm  generates  a solution 
to  (3.1)  for  X = 0,  then  the  original  LCP  is  solved. 

It  is  well-known  (see  [11])  that  this  approach  is  successful  (in  the 
sense  that  it  terminates  either  with  a solution  to  (1.2)  or  with  the 
conclusion  that  the  problem  is  infeasible)  if  the  matrix  L in  the  given 
LCP  is  either  a P-matrix  or  positive  semi -definite . Obviously,  the  matrix  L 
in  the  LCP  (1.2)  need  not  be  a P-matrix  or  positive  semi-definite  in  order 
for  the  LCP  algorithm  to  be  applicable.  In  particular,  if  the  problem  has 
the  property  that  every  pivot  candidate  on  the  diagonal  during  the  application 
of  the  algorithm  is  positive  (this  is  certainly  the  case  when  L is  a 
P-matrix),  then  the  LCP  has  a solution  and  the  algorithm  computes  it. 
Similarly,  if  the  problem  has  the  property  that  whenever  a pivot  candidate 
is  non-positive,  we  can  find  an  "appropriate"  two  by  two  pivot,  or  else  we 
could  declare  the  infeasibility;  then  the  PPP  algorithm  processes  the  LCP 
without  L being  positive  semi-def inite . The  notion  of  the  applicability  of 


Che  PPP  algorithm  may  be  extended  further  by  considerin 
pivots,  but  we  shall  not  consider  it  in  this  paper. 

The  following  conditions  ensure  the  applicability 


Theorem  3.  Consider  a LCP  (1.2)  with  a PQ -matrix  — ^ L.  The  problem  can  be 
solved  by  the  PPP  algorithm  using  solely  the  one  by  one  diagonal  pivots  if 
there  exists  an  m-vector  p satisfying  the  following  conditions: 

(3.2)  r + Xp  > 0 for  some  X > 0 ; 

(3.3)  for  every  index  set  y c {l,...,m}  for  which  L is  nonsingular, 

'iV  = 0 * rkY>  < 0 and  ?kV)  > 0 

hold  for  no  k€{l m}  , where  (y)  denotes  the  principal  pivotal  transform 

with  respect  to  L 

Y Y 


The  proof  of  the  above  result  is  straightforward  and  is  outlined  as  follows: 
Set  up  the  parametric  LCP  (3.1)  using  the  p ; the  condition  (3.2)  ensures 
that  one  may  initiate  the  PPP  algorithm  at  X = X (the  value  of  the  parameter  X 
is  to  be  decreased  to  zero).  The  condition  (3.3)  eliminates  the  possibility 
that  the  k-th  diagonal  element  is  zero  when  it  is  a candidate  for  the  next 
pivot.  The  finiteness  of  the  algorithm  may  be  guaranteed  Dy  usi’r.?  a 
lexicographic  scheme  (see  Cottle  [ 3 ])  or  in  some  cases  by  using  the 
"least-index"  rule  without  lexicography  (see  Chang  [2]  and  Kaneko  [14]). 


5 A real  square  matrix  is  P^  if  all  of  its  principal  minors  are  nonnegative. 


ac  t ..Jcig'' 


We  could  state  a set  of  conditions  ensuring  the  applicability  of 
the  PPP  algorithm  for  a possibly  infeasible  LCP  using  two  by  two  as  well  as 
one  by  one  pivots;  instead  of  discussing  it  in  a fully  general  framework  we 
shall  consider,  in  the  next  section,  a specific  example  of  such  a LCP  arising 
in  a certain  economic  equilibrium  model.  In  the  remainder  of  this  section 
we  shall  identify  LCp's  arising  in  various  applications  satisfying  the 
conditions  in  Theorem  3. 

Example  1 : Consider  the  n by  dn  LCP  (1.1),  where  we  assume  that  (M^ , . . . ,Md) 

has  the  P-property,  and  the  ordinary  LCP  obtained  from  it.  It  was  shown  in 
Kaneko  [12]  that  the  matrix  L is  a p^-matrix  (but  not  a P-matrix)  and  for  p 
such  that  Pj  = 1 , J * 1 n , Pj  ■ 0 , j > n+1 , it  holds  that 

L^^  = 0 implies  that  p£^  * 0,  for  each  y with  det  L^  ^ ^0.  Thus,  the 
conditions  in  Theorem  3 hold.  We  note  also  that  the  PPP  algorithm  applied 
to  the  converted  LCP  can  be  "condensed"  so  that  only  n (rather  than  dn) 
equations  need  be  dealt  with. 

It  follows  from  the  discussion  in  Section  2.2  that  if  the  condition  in 
Theorem  2 is  satisfied,  then  the  solution  of  the  system  of  equations  (2.1)  can 
be  computed  by  using  the  PPP  algorithm  (applied  to  the  corresponding  n by  dn 
LCP).  In  this  connection  we  would  like  to  point  out,  without  elaboration, 
that  the  PPP  algorithm  applied  to  the  n by  dn  LCP  generates  the  same 
solution  path  as  Katzenelson' s algorithm  (see  [15])  applied  to  the  corresponding 
system  of  equations. 


Example  2:  Consider  the  strictly  convex  quadratic  program: 


minimize  f x + £xTFx  subject  to  Ax  < b , x > 0 , 


where  f € Rn  , F is  an  n by  n symmetric,  positive  definite  matrix, 
m ^ _mxii  , , , 

AtR  and  b is  a nonnegative  m-vector.  Notice  that  not  every  feasible, 
strictly  convex  quadratic  program  can  be  put  into  this  form.  It  has  been 
shown  in  Pang  [27]  that  the  LCP  corresponding  to  the  Kuhn-Tucker  conditions 
of  the  above  problem  satisfies  the  conditions  in  Theorem  3. 


Example  3;  The  least  distance  program  (see  Wolfe  [36]  and  Cottle 
and  Djang  ( 5 ])  is  the  problem  of  finding  a point,  with  the  smallest  distance 
from  the  origin,  in  a convex  polyhedron  given  as  a convex  hull  of  m points 
in  Rn.  This  problem  can  be  formulated  as  the  quadratic  program: 

T X X 

minimize  \ x P Px  subject  to  e x = 1 ; x > 0, 


where  P is  an  m by  n matrix  and  e is  the  n-vector  of  one's.  We  note  that 
this  program  is  not  a special  case  of  the  quadratic  program  in  the  previous 
example  (which  requires,  among  other  things,  that  x = 0 be  feasible).  After 
a certain  two  by  two  principal  pivot  performed  in  the  Kuhn-Tucker  condition s, 
the  problem  can  be  transformed  into  a LCP  with  an  n - 1 by  n - 1 matrix 
(see  [ 5]).  By  using  some  special  structures  of  the  problem  and  some  results 
in  [5],  it  can  be  shown  (see  [14])  that  this  LCP  satisfies  the  conditions 
in  Theorem  3 (with  any  positive  p) : in  fact,  it  holds  that  L^^  = 0 implies 

r,^  > 0 for  every  v with  det  L t 0 . 
k = J ' y y 


Example  4;  Consider  the  economic  spatial  equilibrium  problem  with  linear 
prices  mentioned  in  Section  1.  This  can  be  formulated  as  a convex  program 
of  the  form 


(3.4) 


minimize  gxTATDAx  + qTx  subject  to  x > 0 
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where  A is  the  node-arc  incidence  matrix  of  a digraph,  D is  a positive 
diagonal  matrix  and  the  vector  q is  defined  in  terms  of  the  transportation 
costs  and  some  other  constants.  We  note  that  this  program  is  not  a special 
case  of  the  last  two  examples.  Under  the  assumption  that  the  transportation 
costs  satisfy  certain  triangle  inequalities  it  was  proved  in  Pang  and  Lee  [28] 
that  the  LCP  arising  from  this  problem  satisfies  the  conditions: 

L.^  = 0 implies  r,^  > 0 
kk  v k => 

for  every  v with  det  L 0.  It  will  be  shown  in  the  next  section  that 

Y Y 

this  property  carries  over  to  the  case  where  we  allow  the  prices  to  be 
piecewise  linear. 


4. 


4.  A SPECIAL  n BY  dn  LCP 

In  this  section,  we  study  the  following  linear  complementarity  problem: 
given  a positive  integer  -t,  an  n by  n matrix  M;  m by  n matrices  A,N2,...,N^; 

m by  m matrices  B ; an  n-vector  r and  m-vectors  a,  a , ... 

£ or 

with  a >0  for  i=3,...,-L,  find  an  n-vector  x and  m-vectors  y t . . . >y‘t’  such  that 


v i T 1 

(4.1)  u = r + Mx  + r (N  ) y > 0 

i=2 


2 i i 

a + Ax  + Z B y > 0 

i=2 


x > 0 


y2  > 0 


i i 
v = a 


uTx  = (v1)^1  = 0 


y > 0 i = 


i = 


By  decomposing  the  vector  x into  (l-  1)  pieces,  it  is  easy  to  see  that 
this  problem  is  an  (nri-n)  by  (1,  - 1)(m  + n)  LCP  with  the  matrices 


(4.2) 


M1  = / M (N1)1  \ i=2, . . . ,1 


A 3 


Conversely,  an  n by  dn  LCP  of  the  form  (1.1)  is  equivalent  to  a problem 
of  the  form  (4.1). 

The  study  of  the  problem(4.1)  is  motivated  by  its  application  to  the 
piecewise  linear  spatial  equilibrium  model  which  is  of  considerable  importance 
as  briefly  outlined  in  the  Introduction.  In  this  application,  M is  an 
arc-arc  weighted  adjacency  matrix  of  the  (weighted)  digraph  G with  node  set  V 


I 
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denoting  the  various  regions  in  the  market,  arc  set  E denoting  the  routes 
between  regions,  and  (positive)  weights  {b^  ] 3agv  on  the  nodes;  (see  [28]  for 
some  basic  properties  of  such  an  arc-arc  matrix  M) ; each  B1  is  an  identity 


matrix . 


N1  = (D1  - DX)A 


where  is  the  diagonal  matrix  with  positive  diagonal  entries  £b  ^a£\j  » A. 
is  the  node-arc  incidence  matrix  of  the  digraph  G;  r is  an  arc-vector  such 
that  for  each  arc  (cr  , (3)  € E 


(4.3) 


r«(S  ■ a«  - % +c«6 


where  a is  a positive  scalar  and  c is  the  nonnegative  unit  transportation 

“P  l 

cost  from  region  or  to  region  3.  Each  component  of  the  vector  Z y1  where 

i=1 

1 2 2 

y = a - v , represents  the  total  net  import  at  a particular  region.  The 
price  of  the  commodity  at  region  a is  given  by 


p = a - 2 b . y 

« ® i=1  O'1  « 


which  is  a piecewise  linear  function  of  the  total  local  net  import 


i i 

y = Z y by  the  complementarity  conditions  on  the  y 's.  (cf.  (2.5)  and 

“ i=1  “ 

the  discussion  there.)  Each  a^  is  thus  the  equilibrium  price  of  the  commodity 
in  the  absence  of  imports  and  exports. 


I 
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When  -1=1,  the  above  piecewise  linear  model  reduces  to  the  linear  model 
(3,4).  Notice,  however,  that  problem  (4.1)  does  not  in  general  correspond 
to  a convex  quadratic  program. 

Notations . Throughout  the  section,  we  denote  by  M the  matrix 


M (N2)T.  . . . (NV  \ 


A B 


. B 


For  QfC{l,...,n}  and  each  sequence  {3-,..., (3.  } of  disjoint  subsets  of 


{l,...,m}  we  denote  the  corresponding  "principal"  submatrix  of  M 


/ M 


a y 


(n i y ... . (Ng  y ^ 


p2«y 


V 


PP 


. . B 


P P 


by  B(of,3_ , . . . ,3  . ).  Here  3 = U P.  . 
1 * k=2  * 


Definition  1 . We  say  that  the  matrix  M is  consistent  with  M 
in  determinental  signs  if  for  all  a , 3 ^ > • • • >P^>  we  have 

sgn  det  B(or  , 32> . . . ,3  ^ ) = f sgn  det  Ma  a if  a £ 0 

1 otherwise 


where  for  a scalar  x, 


sgn(x)  = 


f’ 


0 

1 


if  x > 0 
if  x = 0 
if  x < 0 


£L . '5*l  _ .*£  ■ - -1. ' 


This  property  of  determinental  consistency  is  motivated  by  the  piecewise 

linear  equilibrium  model  where  it  is  satisfied  easily.  In  general,  if 

2 /_  „ 

(B  , has  the  P-property,  then  the  matrix  M is  consistent  with  M in 

determinental  signs  if  and  only  if  for  each  nonempty  a c {l,...,n} 
sgn  det  (B(a  , P2 , . . . , P^) / (b| p ,...,B  ))  = sgn  det 

2 

where  (B(a  , j^,  • • • (Bg  g ^ denotes  the  Schur  complement 

2 

of  ,#,,,Bq^  ) In  B (of  , ^2  > • • • j 3^)  (see  Co  tt  le  [ 4 1 ) j i • e • 

2 


(B(flf  , P2,-..»P4)/(Bp  p .....B^g  ^)) 


Mof«"((NP2a)  ",,(Np^a)  )(Bp  P 2’*"’BP  P^  Apa 


By  some  easy  calculation,  it  can  be  shown  that  in  the  piecewise  linear 

equilibrium  model,  this  Schur  complement  is  again  an  arc-arc  weighted 

adjacency  matrix  whose  associated  digraph  has  precisely  the  same  structure 

as  that  of  a > only  the  weights  on  some  of  the  nodes  are  suitably  changed. 

According  to  a result  established  in  [28],  the  sign  of  the  determinant  of 

an  arc-arc  matrix  depends  only  on  the  structure  of  the  graph.  Therefore 

the  desired  property  of  determinental  consistency  follows. 

Notice  that  if  M is  consistent  with  M in  determinental  signs,  then 
2 i 

the  matrix  (B  ,...,B  ) has  the  P - property.  The  converse  is  not  always 
true.  However,  we  have 


' 


has  the  p - property  if  and  only  if  M is  a P-matrix  and  M is  consistent 
with  M in  determinental  signs. 


1 Combined  with  Theorem  1 , the  last  proposition  gives  the  following 

j 

; existence  and  uniqueness  theorem  for  problem  (4.1). 

I 

2 3 t 

Theorem  4.  The  problem  (4.1)  has  a unique  solution  for  every  r,a  ,a  , ...,a 

with  a >0  for  i > 3 if  and  only  if  M is  a P-matrix  and  M is  consistent 
with  M in  determinental  signs. 

J 

Remark . Since  the  property  of  determinental  consistency  is  always 
satisfied  in  the  piecewise  linear  equilibrium  model,  the  model  therefore 
j has  a unique  solution  for  all  constant  vectors  if  and  only  if  the  arc-arc 

matrix  is  P,  or  in  other  words,  if  and  only  if  the  network  is  a forest 

(see  [28]). 

' 

As  a consequence  of  Proposition  1 , it  follows  that  if  M is  a P-matrix 
and  if  the  property  of  determinental  consistency  is  satisfied,  then  the  PPP 
algorithm  will  always  compute  the  unique  solution  to  the  problem  (4.1)  by 
performing  only  the  1 x 1 diagonal  pivots.  We  would  like  to  point  out 
nevertheless,  that  in  the  application  of  the  algorithm,  there  is  no  need 
to  introduce  the  matrices  M1.  In  fact,  it  is  easy  to  develop  a condensed 
version  of  the  algorithm  which  operates  on  the  matrix  M.  For  more  details, 
see  [12,  26]. 

In  the  remainder  of  this  section,  we  drop  the  assumption  that  H is  a 
P-matrix  and  study  the  solution  of  (4.1).  To  start,  we  give  an  example 


to  show  that  the  matrix 


obtained  by  writing  (4.1)  in  the  form  (1.2)  need  not  be  positive 
semi-definite  (or  copositive)  even  if  M is  so,  and  if  each  BL  is  the  identity- 
matrix  and  the  property  of  determinental  consistency  is  satisfied. 


With  x 


Hence,  L is  not  copositive,  thus  not  positive  semi -definite 


Proposition  2.  The  matrix  L is  Pn  if  M is  so  and  if  M is  consistent  with  M 


Proof  (for  the  case  K,  = 3) . A typical  principal  submatrix  K of  L has  the 


form  (after  an  appropriate  principal  rearrangement) 
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/M, 


at  a 


<nb  y 


(N2  'T 


3 T 

(N  ) 

Y a 


3 T 

V’  ' 


K = 


S3  a 


M 


3 p3 


P3  P2 


Ms 


M 


Ms 


M2 


M 


p2  e3 


'h  p3 


where  Qf^fl,...,n}  and  y , and  are  mutually  disjoint  subsets  in 

f1,...,m)  . Clearly  det  K is  zero  unless  y = 0 . If  y = 0,  then  by  expanding 

det  K with  respect  to  - IQ  „ , we  have 

p3  p3 


det  K = det  B(cr  , P2  , P3)  >0  . 
by  the  assumptions. 

We  now  establish  that  under  certain  assumptions  which  are  satisfied  in 
the  piecewise  linear  equilibrium  model,  the  PPP  algorithm  will  successfully 
process  the  problem  (4.1). 

Theorem  5.  Suppose  that  the  conditions  below  are  satisfied: 

(4,5)  M * eVe  where  E 6 Rm  X n and  D1  € Rm  X m with  D1  being  positive 
definite; 
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(4.6)  A = CE  and  N1  = D1 E for  i = 2 4 where  C 6RmXm  and  D1  GRmXm  ; 

(4.7)  The  matrix  M is  consistent  with  M in  determinental  signs. 

5/  T T T 

Then  by  choosing  the  parametric  vector—  to  be  (en  , , 0,...,0)  , the  PPP 

algorithm  will  in  a finite  number  of  pivots,  either  terminate  with  a 
complementary  solution  to  the  problem  (4.1),  or  indicate  that  the  problem 
has  no  feasible  solution. 

Before  proving  this  theorem,  we  point  out  several  remarks.  First, 
assumption  (4.5)  implies  that  M is  positive  semi-definite  but  not  necessarily 
symmetric  or  positive  definite.  Second,  as  Example  5 above  shows, 
these  assumptions  do  not  imply  that  the  matrix  L is  positive  semi -definite 
or  copositive.  Finally,  except  for  the  assumed  property  on  D*,  the  matrices  DJ , 
C and  E are  otherwise  arbitrary. 

We  need  the  lemma  below  to  prove  the  theorem.  It  follows  from 
Property  2.3  in  [12]. 

Lemma  1 . Under  the  assumptions  of  Theorem  5,  the  following  assertion  holds: 

If  a E >n}  fs  nonempty  such  that  M^a  is  nonsingular,  and  if  {°2> ,{3^} 

is  a family  of  mutually  disjoint  subsets  of  {l,...,m},  then  for  2 < i < j < 4, 
the  matrix 


® (®  >32 » • • • > B^)  ® (®  ,82 9 j 

4 

is  P.  Here,  for  p » U 9.  , 

k-2  * 


•»4> 


5/  By  e^  we  mean  the  k-vector  of  ones. 


® (**  9 ^2  5 • • • * > • • • > P j > • • • > P ^ = 


2 T 

<V>  " 


••  <*)T 


<»>T- 

J 


...  Bi 


• (N3  a)T  ■ 


...  B" 


Proof  of  Theorem  5.  To  simplify  the  notations,  we  shall  prove  the 
theorem  for  1=3.  The  general  case  can  be  treated  in  a similar  way. 
Let  ac  n)  , &2  and  P3  E be  index  sets  such  that 


(i)  P2  n p3  = 0 


x » yf  HQ  and  are  currently  basic. 
® °2UP3  P-* 


The  fact  that  the  set  of  basic  y -variables  always  contains  that  of 
3 

basic  y -variables  may  be  justified  easily.  (See  [12,  26 ] e.g.)  The 
current  canonical  system  has  the  form 
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Constant  Parametric  x^ 

Column  Column 
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M M . 

a or  a 6 


i i 


6 

3 

'*3 
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M,  M 
6 a 6 6 


a -2 

Aa  . - Bn  B 

0 3®  ® 3®  3 
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B2 
Y Y 
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Consequently,  the  next  diagonal  pivot  entry  will  be  positive  unless  it  occurs 
at  a Ug-row.  Suppose  now  that  the  next  diagonal  pivot  entry  does  occur  at 

A 

a u. -row  with  i £ a and  m. . is  zero.  Then  it  follows  that 
i 11 


det  B(of'  , &2  , p3)  =0 


where  a'  = orU  {i } . By  assumption  (4.7),  H#I  , is  singular.  Therefore 
E.^,  has  linearly  dependent  columns.  This  implies  that  for  some  vector 


f 0 , we  have 
or  ’ 


"M-  JWIW,IUJ-1 11  * 
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and 


M = M 
6 i 6i 


- (M  (N3  )T  (N3  )T 
^ 6a  '•w336;  ^NB26; 


, M . , 
/ er  i \ 


031  i 

U ,/ 


= 0 


Consequently,  in  the  current  canonical  system,  except  for  those  x -rows 

0 1 9 

the  entire  xfi -column  is  zero.  Furthermore,  we  have 

fli«  “ <Ml«WB33i)T  (Mo1)T 

Therefore,  one  can  perform  a 2x2  block  pivot  unless  f >0.  If  f >0 

Of  = 0/  = 9 

then  the  homogeneous  problem 


l'\ 

-/ 

2 

V 

\vV 

\ 

/ N (U2)T  <N3)T 

2 _3 


B 


I *\i° 


B 

0 


y 

\ 3 
\ y 


r\i0 

2 I 

y 

W3/ 


T . 2.T  2 . 3.T  3 

u x = (v  ) y = (v  ) y =0 


2 3 

has  a nonzero  complementary  solution  given  by  y = y =0  and 


Xy  ■ ftf  » Xj  - 1 and  x^  = 0 for  j £ a'  . 


^ ij  ^ 2 m o 

Now  x Mx  = 0 implies  N x ■ N x=  0 so  that  if  X*  is  the  current  value  of  the 


parameter,  we  have 


* 
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0 = xTu*  = x T (r  + Xe)  + xT  [Mx*  + E(N1)T(y1)*] 

~ T T ~ — T 

» x r + \*e  x > x r 

because  X*  is  positive.  Consequently  the  system 

u - r + Mx  + S(Nl)V  >0  , x > 0 y1  > 0 

and  therefore  the  problem  (4.1),  is  infeasible.  This  completes  the  proof 
of  the  theorem. 

The  last  part  of  the  proof  also  establishes 

Corollary  1.  Under  the  assumptions  of  Theorem  5,  if  either  one  of  the 
two  systems 

(i)  • u = r + Mx  > 0 , x > 0 

(ii)  u = r + Mx  + Z(Ni)Tyi  >0  x > 0 , yL>0 

is  consistent,  then  the  PPP  algorithm  will  always  compute  a solution  to 
problem  (4.1). 

Remark.  Theorem  4.5  remains  valid  if  M is  given  by 
M = eVe  + Dl+  *E 

provided  that  the  product  D^+^E  is  positive  semi -definite. 

In  general,  the  conditions  in  Theorem  5 are  not  sufficient  to 

guarantee  that  the  PPP  will  terminate  with  a solution  to  the  problem  (4.1). 
In  what  follows  we  present  a sufficient  condition  which  together  with  those 
of  Theorem  5,  will  imply  that  the  conditions  in  Theorem  3 are  satisfied.  In 
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r : 


j 

1 

. 

! 

I 

1 

I 


i 


particular,  the  problem  (4.1)  will  then  have  a solution  which  can  be  computed 
by  the  PPP  algorithm  by  performing  only  the  1 x 1 diagonal  pivots. 


Theorem  6.  Suppose  that  the  assumptions  (4.5)  - (4.7)  in  Theorem  5 are 
satisfied.  Futhermore,  suppose  that  there  exists  an  n-vector  s satisfying 


(4.9)  r + \s  > 0 for  some  scalar  X > 0 

(4.10)  For  every  arc  [l,...,n}and  ilfar  such  that  the  columns  in  E.  are 
linearly  independent,  it  holds  for  no  vector  f such  that 

E.,  = E.  f r.  - rTf  < 0 and  s.  - sTf  > 0 . 
i Cfor  larar  iaor 

Then  conditions  (3.2)  and  (3.3)  in  Theorem  3 are  satisfied  for  the  problem 
(4.1)  when  it  is  transformed  as  an  ordinary  LCP  of  the  form  (1.2)  with  the 
matrix  L given  by  (4.4). 


Proof.  Observe  that  in  the  current  canonical  system  (given  in  the  proof  of 
Theorem  5),  the  only  place  where  condition  (3.3)  could  possibly  be  violated 
is  in  a u^-row  with  ijEor  and  = 0.  By  the  argument  used  in  the  proof  of 
Theorem  5,  it  is  not  difficult  to  show  that  the  current  values  of  the 
corresponding  r-  and  s -components  are  given  by 


r - rTf 
i a or 


and 


s , = s - s f 
i i a a 


respectively.  Here  f is  given  in  (4.8).  From  the  given  assumption  (4.10), 
the  desired  conclusion  follows  readily. 

As  a concluding  remark,  we  point  that  in  the  piecewise  linear 
equilibrium  model,  the  condition 


E. 


i 


= E 


. f 
a a 


implies  that  the  arcs  indexed  by  arl)  { i } form  a minimal  cycle.  As  proved 
in  1 28] , this  cycling  condition  in  turn  implies  that  the  current  i-component 
of  the  r-vector  (i.e.,  r^)  is  nonnegative,  provided  that  the  transportation 

costs  c (cf.  (4.3))  are  nonnegative  and  satisfy  the  triangle  inequality 

a b 


c _ + c„  < c 
<*9  9 Y = Of  Y 


Consequently,  under  these  conditions  on  the  transportation  costs,  the  last 
theorem  implies  that  the  PPP  algorithm  can  be  successfully  used  to  compute 
a solution  to  the  piecewise  linear  equilibrium  model. 
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